Prom ageing to immortality: cluster growth in stirred colloidal solutions 
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Our model describes cluster aggregation in a stirred colloidal solution Interacting clusters com- 
pete for growth in this 'winner-takes-all' model; for finite assemblies, the largest cluster always 
wins, i.e. there is a uniform sediment. In mean-field, the model exhibits glassy dynamics, with two 
well-separated time scales, corresponding to individual and collective behaviour; the survival prob- 
ability of a cluster eventually falls off according to a universal law (int)" 1 ^ 2 . In finite dimensions, 
the glassiness is enhanced: the dynamics manifests both ageing and met ast ability, where pattern 
formation is manifested in each metastable state by a fraction of immortal clusters. 

PACS numbers: 68.43.Jk, 47.54.+r, 89.75.-k 



Systems that are far from equilibrium exhibit remark- 
able physics: some commonly known examples in the 
context of glassy systems [lj are ageing and metasta- 
bility. There are, however, much simpler consequences 
of nonequilibrium dynamics; in this Letter, we present a 
model of cluster growth which is inherently disequilibrat- 
ing. Here, instead of available masses being equally dis- 
tributed between competing clusters, the largest cluster 
always wins. Somewhat surprisingly,the two-step relax- 
ation and ageing characteristic of glassy dynamics is man- 
ifested for infinite systems. Additionally, in finite dimen- 
sions, the system relaxes asymptotically to a metastable 
state in a complex energy landscape; in each such state, a 
fraction of immortal clusters from nontrivial spatial pat- 
terns. 

In our model, n pointlike, immobile clusters with time- 
dependent masses rrii(t) for i = 1, . . . , n, evolve according 
to the following equations: 
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These dynamical equations were originally written to 
model the kinetics of black hole growth in a radiation 
fluid 0, E»]- The present model is, however, of poten- 
tial interest in many other situations; its dynamical be- 
haviour in fact illustrates the rich-get-richer principle 
in economics Q, and its form bears some resemblance 
to Lotka-Volterra type predator-prey models || in bio- 
physics. In this Letter, we visualise Eqs. represent- 
ing a simple model of aggregating clusters in a stirred 
colloidal system 6]. In our toy model, clusters grow by 
accreting mass from the fluid and from other clusters; 
they can also dissolve away if they are small, as the fluid 
is stirred. The gain term in Q is the sum of the free 
rate for an isolated cluster, proportional to the parame- 
ter a > 1/2, and of constributions from all other clusters 
via the fluid, with a coupling gij between clusters i and 
j; the loss term is taken to be inversely proportional to 
the cluster mass. The stirring of the fluid results, in this 
simple-minded picture, in the explicit time dependences 



in 0. 

We first recall 0, the- one-cluster result. (For con- 
venience, we work with reduced time s = In ^ (to renor- 
malise away the effect of initial time to), reduced masses 

2 

x i — jt?2 and square masses y, = x 2 = -j*-.) Large clus- 
ters, whose initial mass yo is greater than some threshold 
y*, with y*(to) = ( 2^-1 ) ' are survivors: they survive 
and keep on growing forever. Small clusters, whose ini- 
tial mass is below this threshold dissolve and die out in 
a finite time. In a very dilute sediment, we would thus 
expect globules of matter to be suspended indefinitely, 
and to grow forever. 

As the concentration of the solute is increased, inter- 
actions between the colloidal clusters become significant. 
We consider first two interacting clusters with interac- 
tion strength g. If their initial masses are exactly equal, 
this equality is maintained by symmetry forever. The 
reduced mass x(s) of each mass then obeys: 



(2a - l)x 2 - 2 - gx 3 
2x(l + gx) 



(2) 



The fixed points of|21are given by (2a — l)x 2 — 2— gx 3 = 0; 

there is a critical value g c = {^—^f—^j which sepa- 
rates two kinds of behaviour. For large couplings g > g c 
, there is no fixed point; physically this implies that the 
clusters feed on each other and disappear quickly. For 
small couplings g < g c , there are two positive fixed points 



yV 2 < sc(i) (unstable) < (iy^) 1 ^ 2 < x^) (stable). (3) 

Small clusters, such that their mass Xo < Xm, are dy- 
namically attracted by a; — 0; thus lighter clusters dis- 
solve away rapidly. Large clusters, such that their mass 
xq > X(i), are dynamically attracted by X(2y. heav- 
ier clusters thus grow forever according to the law as 
m(t) X(2)t x l 2 . Here the interaction is small enough 
that symbiosis is achieved; the clusters feed on each other 
so as to increase, rather than deplete, mutual growth. 
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For two clusters with initially unequal masses, any 
small initial mass difference always diverges exponen- 
tially at early times: later stages of the dynamics cannot 
be described in closed form. The details of the transient 
behaviour can be found in 0. However, the asymptotic 
behaviour is such that the largest cluster wins: survival 
of the biggest is therefore the single generic scenario for 
two clusters with unequal initial masses. The surviving 
large cluster may then cither also disappear in a finite 
time, or survive and grow forever, depending on whether 
its mass at the time of the lighter cluster's death is less 
or greater than the threshold y*. This generalises easily 
to any finite number n > 2 of colloidal clusters; in every 
case, our model predits that there will at most be one 
large sediment formed 0. 

To explore an infinite assembly of interacting col- 
loidal clusters, we explore the mean- field behaviour of 
the model; each cluster is connected to every other by a 
weak interaction g = — . When n — > oo limit at fixed g, 
the socalled thermodynamic limit, one obtains |2J: 



y '(s) = 7(s)y(s) 



(4) 



for the reduced square mass y(s) of any of the clusters. 
Eq. 10} can be solved only formally, as it is self-consistent 
0. 

Things simplify considerably when the rescaled cou- 
pling 3 is small; remarkably, a glassy dynamics 00 with 
two-step relaxation is observed. In Stage I, the clusters 
behave as if they were isolated: this corresponds to indi- 
vidual behaviour. The survivors of this stage are clusters 
whose initial masses exceed the threshold y*, exactly as 
in the one-cluster case recalled above. In Stage II, all 
the survivors interact with each other; the dynamics is 
thus collective, and turns out to be slow 0. All but the 
largest cluster eventually die out during this stage, thus 
in terms of experimental predictions , again predicting 
a uniform asymptotic sediment. This weakly interact- 
ing mean-field regime of our model of interacting col- 
loids clusters also exhibits other characteristic features 
of glassy systems such as ageing; this originates in 
the presence of two well-separated time scales of fast and 
slow dynamics, whose ratio grows as 1/g 2 [3. 

For an exponential distribution of initial masses, in the 
late times of Stage II, the survival probability decays as 



S(t) 



with 



H^fcini 
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irrespective of a, g and the parameters of the exponen- 
tial distribution. Also, the mean mass of the surviving 
clusters grows as 

1/2 
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The universality inherent in the scaling results ©-0 
is unusual, because it includes the prefactor C, which is 
itself independent of the details of the initial distribution 
P(yo) of square masses. It can in fact be shown that C 
only depends on the tail exponent of this distribution in 
the vicinity of its upper bound y ma x, i-e. whether the 
initial distribution of masses is bounded or not. The in- 
terested reader is referred to for further details; here 
we simply point out that this striking universality is a 
major result of our present work. Also, the logarith- 
mic behaviour seen is yet another signature of glassy 
dynamics. Returning to the physical system, this pre- 
dicts that in a very weakly interacting system of colloidal 
particles, there is a slow freezing in of fluid disorder |||, 
resulting in a uniform sediment at asymptotically long 
times. 

Finally, we study a lattice version of the model ( we 
choose the chain (D — 1), the square lattice (D = 2), 
and the cubic lattice (D = 3)) in order to probe the 
effect of fluctuations. Clusters now sit on the vertices 
n of a regular lattice, with nearest-neighbour interaction 
g. In the limit of weak coupling (<? <C 1), our dynamical 
equations are: 



2a - 1 
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where m runs over the z = 2D nearest neighbours of site 
n, 

The dynamics generated by (jHJ , again consists of two 
successive well-separated stages. Fast individual dynam- 
ics are exhibited in Stage I, where the mass of each clus- 
ter evolves as if it were isolated. As in the mean-field 
case, the survival probability S(s) decays rather fast from 
5(0) = 1 to a plateau value The effects of going 

beyond mean field are only palpable in Stage II, where 
interactions become relevant and lead to a slow dynam- 
ics which is now very different from the mean-field sce- 
nario. The survival probability S(s) in fact decays from 
its plateau value Sm to a non-trivial limiting value S(oo)j 
reflective of metastability. The effect of fluctuations is al- 
ready palpable: unlike the mean field result, this predicts 
that a, finite number of clusters will survive. We elaborate 
on this below. 

Consider © for two neighbouring clusters n and m 
which have both survived Stage I. The contribution of 
cluster m to the large parenthesis in the right-hand side 
of ||SJ| is proportional to agx m . In the absence of cou- 
pling, we have x m ~ e^ 2a_1 ^ s / 2 00. The characteris- 
tic time scale of Stage II, that is the time at which in- 
teractions become significant, is reached when the prod- 
uct gx m becomes of order unity: 



2a 



1 9 



(9) 



(7) 



i.e. t c ^tog 2 /( 2q !). Thus, in this weak-coupling limit, 
the separation of time scales between fast individual and 
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slow collective dynamics is parametrically large. Figure^ 
illustrates this two-step relaxation in the decay of the sur- 
vival probability 5(s) in one dimension. Both stages of 
the dynamics appear clearly on the plot; different values 
of the interaction cause the system to age differently. In 
each plateau is reached at Sm — 0.8; however, 

the weaker the interaction, the longer the system takes 
to reach the asymptotic state, which occurs at a non- 
trivial limit survival probability S(oo) ~ 0.4134. Since 
each curve corresponds to a decade shift in interaction 
strength g, it is shifted in terms of the onset time s c of 
slow dynamics by 2 In 10 (thick bar), in excellent agree- 
ment with the estimate (j^J. 
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FIG. 1: Plot of the survival probability S(s) on the chain 
with S m = 0.8. Left to right: Full line: g = 10" 3 . Dashed 
line: g = 10 -4 . Long-dashed line: g = 1CP 5 . Dash-dotted 
line: g = 10" 6 . The thick bar has length 2 In 10 = 4.605 (see 
text). 

At the end of Stage II, the system is left in a non- 
trivial attractor, which consists of a pattern where each 
cluster is isolated: it is therefore a survivor and keeps 
growing forever. We call these attractors metastable 
states, since they form valleys in the existing random en- 
ergy landscape; the particular metastable state chosen by 
the system (corresponding to a particular choice of pat- 
tern) is the one which can most easily be reached in this 
landscapeP,0,0. The number Af of these states gener- 
ically grows exponentially with the system size (number 
of sites) N: 

A/"~exp(AT£). (io) 

with £ the configurational entropy or complexity. 

The limit survival probability 5(oo) (Figure Q is just 
the density of a typical attractor, i.e., the fraction of 
the initial clusters which survive forever. It obeys the 
inequalities 

S(oo) < 5(1), 5(oo) < 1/2. (11) 

The first inequality expresses the fact that clusters gener- 
ically disappear: the difference 1 — 5(1) (resp. S^—S^^) 
is the fraction of clusters which are dissolved out during 



Stage I (resp. Stage II). The second inequality is a conse- 
quence of the fact that each surviving cluster is isolated: 
the densest configuration for which this is the case is 
when either of the two sublattices is occupied, at which 
point the density is exactly 1/2. This value 1/2 of the 
highest density holds for the large family of so-called bi- 
partite lattices, which includes the hypercubic lattices we 
have considered here (but does not, for example, include 
the triangular lattice). 

For a given class of initial mass distributions, the limit 
survival probability St^) is a monotonically increasing 
function of the plateau value S'(i); the more the number 
of survivors after Stage I, the higher will evidently be 
the number of immortal clusters. For SVqo) = 0, Sm 
is trivially 0; as the 5m — + 1 limit is approached, the 
non-trivial maximum value 5( oc ) max < 1/2 is reached. 
Additionally, when 5m is small, it can be shown that 
^(oo) is also small, and that it depends on 5m alone. 
This is shown below. 

We define a supercluster as a set of k > 1 connected 
clusters which have survived Stage I, such that all their 
neighbours have disappeared during Stage I. The fate of 
superclusters depends on their size k as follows. 

★ k = 1: If a supercluster consists of a single iso- 
lated cluster, it is a survivor, because its mass ex- 
ceeds the survival threshold y*. For z — 2D and 
independently of initial mass distributions, a su- 
percluster with k = 1 occurs with density p\ = 
5(1) (1 — 5(i)) 2£) . This corresponds to the survival 
of one cluster after, and the disappearance of its 
2D neighbours during, Stage I. 

★ k = 2: If a supercluster consists of a pair of neigh- 
bouring clusters (represented as ••) both clusters 
evolve according to the dynamics described above: 
the smaller dies out, while the larger is a survivor. 
We are thus left with »o or o» in the late stages 
of the dynamics. Such an event takes place with 
density p 2 = 5f 1) (l-5(i ) ) 2 ( 21? - 1 ). 

★ k > 3: If three or more surviving clusters form 
a supercluster, they may a priori have more than 
one possible fate. Consider for instance a linear 
supercluster of three clusters (•••)■ If the middle 
one disappears first (»o»), the two end ones are iso- 
lated, and both will be survivors. If one of the end 
ones disappears first (e.g. ••°), the other two form 
an interacting pair, and only the larger of those two 
will survive forever (e.g. too). The pattern of the 
survivors, and even their number, therefore cannot 
be predicted a priori. 

The above enumeration implies 5(00) = Pi + P2/2 + • • • , 
where the dots stand for the unknown contribution of 
superclusters with k > 3. As p\ ~ 5(i), P2 ~ <Sm, an d 
so on, we are left with the expansion 

S(oo) = 5(i) — D 5( X ) H (12) 
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The dependence of SVoo) on details of the initial mass 
distribution at fixed 5(i) therefore only appears at order 
Shy These results apply to the dilute limit, when few 
clusters survive stage I. 

In the limit that there are many clusters which sur- 
vive the fast dynamics of stage I, i.e. Sen — > 1, the 
limit survival probability, as mentioned above, reaches 
a non-trivial maximum value iSVooWax < 1/2- This de- 
pends very weakly on the mass distribution; for instance, 



in one dimension one has Si 



ncntial distribution and S, 



(oo)i] 



(oo)n 



« 0.441 for an expo- 
0.446 for a uniform 



distribution. We present below a way of visualising these 
dense distributions of immortal clusters. 

If (S^) = 1/2) on, say, a square lattice, (i.e. the high- 
est density of immortal clusters is reached) , there are only 
two possible 'ground-state' configurations of the system. 
These correspond to the full occupation of one of the two 
sublattices, with its counterpart being completely empty. 
In this limit, the two possible patterns of immortal clus- 
ters are each perfect checkerboards of one of two possible 
parities. 

At every site n with integer co-ordinates 
(ni, ri2, • • ■ , no) we define the survival index 



1 if the cluster at site n is a survivor 
= else (13) 

and the checkerboard index 



0n = (-l) 



fn+niH Vnu 



(14) 



The survival index depicts very simply the pattern of sur- 
viving clusters surrounded by empty sites: The checker- 
board index, on the other hand, represents, for each 
site, the local choice of one of the two symmetry-related 
'ground states', i.e., of one of the two sublattices. This is 
easiest to understand using a one-dimensional example: 

the two ground states are H 1 1 or — | 1 — ••• 

All the <j) n are equal to —1 in the first pattern, and equal 
to +1 in the second pattern. The checkerboard index 
4> n thus classifies each site according to the parity of the 
particular ground state selected locally at this site. 

For generic initial conditions, i.e. a random distribu- 
tion of initial masses, the immortal sites will evidently 
not form a perfect checkerboard. However, if the initial 
masses are large enough, the number of survivors of Stage 
I dynamics will be large, and the corresponding survival 
probability Sm after Stage I close to unity. In this limit, 
the asymptotic survival probability SVqq) will be close to 
its 'ideal' value of 1/2. The resulting pattern will now 
exhibit a local checkerboard structure, with frozen-in de- 
fects between patterns of different parities; the random 
structure of defects is of course entirely inherited from the 
(random) initial mass distribution, since the dynamics is 
fully deterministic. This is evident from Figure [21 which 
shows a map of the survival index and of the checker- 
board index for the same attractor for a particular sample 
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FIG. 2: Two complementary representations of a typical pat- 
tern of surviving clusters on a 40 2 sample of the square lat- 
tice, with 5(i) = 0.9, so that SVoo) ~ 0.371. Left: Map of 
the survival index. Black squares represent immortal sites for 
which a n — 1, while white squares represent dead sites for 
which (j n = 0. Right: Map of the checkerboard index. Black 
squares represent n = +1, while white squares represent 
(t> n = -1. 



of the square lattice. The immortal (black) clusters in the 
left-hand part of the figure are surrounded by rivulets of 
voids, which are a consequence of initial conditions; in the 
right-hand figure, the deviation from a perfect checker- 
board structure (all black or all white) is made clearer. 
The patterns make it clear that neighbouring sites must 
be fully anticorrelatedm because each immortal cluster 
is surrounded by voids. However, at least close to the 
limit /SVoo) = 1/2, immortal sites are very likely to have 
next-nearest neighbours which are likewise immortal and 
massive. The detailed examination of survival and mass 
correlation functions made in a longer paper confirms 
these expectations. 

To conclude, we have presented in this Letter a very 
simple 'winner-takes-all' model of cluster aggregation in 
stirred colloidal systems. Both mean-field and finite di- 
mensional explorations of this model show a striking and 
a priori unexpected glassy phase; the system of inter- 
acting clusters shows ageing until it reaches its asymp- 
totic state. The inclusion of fluctuations in the model via 
a finite-dimensional approach causes the replacement of 
the somewhat staid mean-field behaviour (which predicts 
a uniform sediment composed of one cluster at most) by 
something far more exciting; a random-energy landscape 
emerges, with many possible minima as its metastable 
states, and the system descends to the most accessible 
one. Each such metastable state is a complex pattern of 
isolated clusters, each of which, by virtue of its isolation, 
is immortal. 
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